Fictitious flux confinement: magnetic pairing in coupled spin chains or planes 
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The spinon and holon excitations of two-leg Heisenberg or lightly doped t-J ladders are shown to be 
bound in pairs by string confinement forces given approximately by the antiferromagnetic exchange 
energy across the rungs, F = J±{SiSj)±/b. These forces originate from the fictitious flux tubes 
associated with the half-fermi statistics of the excitations. It is conjectured that similar confinement 
forces, determined by the antiferromagnetic exchange energy across the layers, are responsible for 
the spin gap and the pairing of charge carriers in CuO superconductors. 

PACS numbers: 71.10.-w, 75.10.-b, 74.72.-h, 74.20.-z 

1. Two-leg Heisenberg or lightly doped t-J ladders j|] 
can, according to a recent proposal (2| , be described ap- 
proximately by an RVB type spin liquid. To be precise, 
one assumes a magnetic tight-binding model on the lad- 
der, with flux TT per plaquet and hopping magnitudes 
t along the chains and t± across the rungs. One then 
fills the lower band twice, once with up-spin and once 
with down-spin electrons. Upon elimination of doubly 
occupied sites via Gutzwiller projection, one obtains a 
reasonable approximation to the ground state of the 
Heisenberg ladder with J±/ J = i±/i. 

Spinon and holon excitations for this liquid may be 
created either via Anderson's projection technique ^ or 
via midgap states [|| . As the topology of the ladder dic- 
tates that midgap states can only be created in pairs, the 
second method automatically yields spinon-spinon bound 
states (magnons) rather then isolated spinous, which re- 
flects the fact that the spinous or holons are confined in 
pairs. 

In the first part of this letter, I will explain the form 
and origin of the confinement forces between the spinous 
and holons of the ladder, and calculate the spinon mass 
and the bound state resonances based on a heuristic iden- 
tification of the spin gap of the Heisenberg ladder as the 
zero point energy of the string oscillator. In the second 
part, I will postulate similar confinement forces in sys- 
tems of (weakly) coupled magnetic planes, make some 
assumptions regarding both the nature of the spin liq- 
uid in the planes and the confinement forces due to the 
interplane coupling, and obtain an estimate for the spin 
and charge gaps in bilayer CuO superconductors (25 meV 
and 28 meV, respectively, for YBa2Cu306+i;). I conclude 
with several comments on high- Tc superconductivity. 

2. In order to determine the functional form and 
strength of the confining potential between the spinon 
or holon excitations of the two-leg Heisenberg ladder, we 
create two holons at sites i and j via Anderson's projec- 
tion technique H , 
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Number of rungs between the two holons 

FIG. 1. The confinement energy, defined here as the energy 
of a pair of holons measured relative to the exact ground state 
energy of a pair of holes, as a function of the number of rungs 
between the holons or holes for a 2 x 8 ladder with open 
boundary conditions and i±/i — 0.75(Jj_/J)°'® for J±/J < 1 
and i±/i — {J±/jy^'^ for J±/J > 2, i.e. such that the ratios 
{SiSj)±/{SiSj}i\ of the trial wave functions match those of 
the exact ground states. One holon was localized at the end 
of the ladder and the other localized at various positions on 
the same chain, as illustrated in fig. 3a; a finite size correction 
has been taken into account for J±/J = 0.2 and 0.1. 
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where IV'Sd) is the Slater determinant ground state ob- 
tained by filling the lower magnetic tight-binding band 
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FIG. 2. The string tension in units of J±_ given by the slope 
of the dotted lines fitted though the data points in fig. 1 in 
comparison with the spin correlations |(5'iS'-,)i| across the 
rungs in the ground state as taken from table I of [2] for dif- 
ferent ratios Jx/J. The discrepancy around Jx/J = 2 arises 
from the enhanced correlations KSiS^)||| along the chains in 
the presence of the holons, which have been neglected in (2). 
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FIG. 3. Spin correlations (SiSj) between nearest neighbors 
in the presence of two holons compared with two holes for 
two representative configurations of a 2 x 8 ladder with open 
boundary conditions and (a) J±/J = 0.5 or (b) J±/J = 5. 
The conventions are adopted from White and Scalapino [6] 
the thickness of the lines is proportional to KSiSj)!; solid 
lines indicate antiferromagnetic, dotted lines ferromagnetic 
correlations. The state containing two holons was obtained 
from Gutzwiller projected magnetic bands using Anderson's 
method for constructing spinous; the state containing two 
holes is just the exact ground state of the Heisenberg ladder 
with two static vacancies at the positions indicated. 



twice, and numerically compare the energy expectation 
value of this configuration to the energy of the exact 
ground state for a Heisenberg ladder with two station- 
ary holes at these sites The results are shown fig. |l| 
and |; we find a linear potential 
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where x is the distance between the spinous or holons in 
the direction of the chains, (8^8^)^ the spin correlation 
across the rungs, and b the bondlength along the chains. 
The confinement energy is hence approximately equal to 
the antiferromagnetic exchange energy across the rungs 
between the spinous; the reason for this emerges from 
a comparison of the individual spin correlations on each 
link, as reproduced in fig. ^ for two typical configura- 
tions: an invisible string between the holons destroys the 
antiferromagnetic correlations on all the rungs between 
them. 

The origin of this string is explained in fig. ^: the fic- 
titious flux tube associated with the half-fermi statistics 
of the spinous or holons , which manifests itself in 
an adjustment by tt |^J^ of the fictitious fiux through the 
adjacent plaquets in the magnetic tight-binding model 
before Gutzwiller projection, effectively annihilates the 
hopping terms on the rungs between them. 

3. If we were to know the effective mass of the spinon, 
we could calculate the resonances of the string oscillator 
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Conventions for hopping phases across links: 
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FIG. 4. Magnetic tight-binding configurations before 
Gutzwiller projection for (a) the ground state of the t-J lad- 
der at half-filling and (b,c) in the presence of two holons at a 
distance of 4 lattice spacings. The holons in (b) were created 
following the procedure for the Kalmeyer-Laughlin chiral spin 
liquid [7] suggested in [2] after mapping the flux ladder into 
a flux lattice subject to a periodic boundary condition with 
a periodicity of only two lattice spacings in y-direction. The 
lattice is subsequently reconverted into a ladder (c). Note 
that the Dirac string annihilates all the hopping terms across 
the rungs between the holons, while the hopping magnitudes 
along the chains remain unaffected. 
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where mrod = ^"^sp is the reduced mass of a spinon pair, 
and compare the ground state energy to the spin gap in 
the Heisenberg ladder, which is just the energy required 
to create a spinon-spinon bound state; apart from a pos- 
sible correction due to a chemical potential for spinous 
which we neglect, these should be equal. In the present 
case, we use the know value A « J±/2 for the spin gap 
pO| to calculate the spinon mass and the spinon-spinon 
resonances in the weak coupling regime J± < J. In one 
dimension, the solutions of are given in terms of the 
Airy function Ai(a;) 

M^) = Ai (^M _ A„^ , Er, = FxoA„, (4) 

where xq — l/(2mrod-F')^/'^ is the characteristic length- 
scale of the oscillator, and A„ are the extrema or zeros of 
Ai(— x) for n even or n odd, respectively, which are listed 
in table ||. Equating Eq = A, we obtain 
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where we have approximated |(SiSj}_L| ~ 0.62J±/J ac- 
cording to table I of dl for weakly coupled ladders. The 
spacing of the eigenvalues A„ implies that the internal 
resonance frequencies of the spinon-spinon bound states 
or magnons are comparable but slightly higher than the 
energies required to create a second or third magnon; 
this proximity and the fact that they presumably decay 
rapidly into several magnons renders a potential obser- 
vation of the resonances difficult if not impossible. 

It is now easy to estimate the size of the magnon. Using 
Ai(|a;| — Aq) w Ai(— Aq) exp(— ia;^), we write the ground 
state 

Mx) = exp(-— ^) with ^ = J^xq = 0.97—5. (6) 
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This result illustrates why the spin gap in the weak cou- 
pling regime can be ^J± while the antiferromagnetic ex- 
change energy across each rung is only of order ^J'^/J: 
the number of decorrelated rungs is of order J/J±. 

This calculation can also be applied to spinon-holon 
bound states (holes) in the ladder; the only difference is 
that the reduced mass in and (^) is replaced by 

1 1 1 i„ 1 
= \ , where Wh = ^, (7) 

is the effective mass of the holon; according to table II 
of [§, teff = 5£^t|| = 0.77t for = J and t^g = 0.95 1 
in the weak coupling limit J± <C J. The values for the 
energies En and the size ^ of the bound state are those 
given in and above multiplied by 
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TABLE I. The lowest (dimensionless) energy eigenvalues 
for the linear oscillator in one and two dimensions. 
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4. I will now turn to the more speculative part of this 
letter, and explain part of my thinking on CuO super- 
conductivity. To begin with, I make the following as- 
sumptions: 1. The ground state of the two-dimensional 
t-J model at the relevant hole dopings is a spin liquid, 
which supports spinous and holons as elementary exci- 
tations 1^. (My understanding jlj] of this liquid is that 
it is a liquid in both the spin degrees of freedom and 
in the non-relativistic plaquet chiralities ; this chiral- 
ity liquid may be seen as a significant generalization of 
Laughlin's chiral spin liquid which is a liquid in the 
spins but effectively aligns the chiralities and thereby vi- 
olates the discrete symmetries P and T. These symme- 
tries are preserved in my construction. The spinon and 
holon excitations supported by the chirality liquid carry 
a chirality quantum number, which can be + or — ; this 
number determines the sign of the winding phases asso- 
ciated with their half-fcrmi statistics.) 2. The spinon 
and holon masses in a system of coupled t-J planes are 
comparable to their values in a system of coupled chains 
estimated above. 3. The analysis above, including all 
the implicit assumptions made, is correct. 

The main difference between fictitious flux confinement 
in a system of coupled planes as compared to coupled 
chains is that the one-dimensional array of decorrelated 
rungs is replaced by a puddle of decorrelated interplane 
links. To see this, imagine several ladders (which are 
not necessarily straight) embedded into a system of cou- 
pled planes such that the rungs align with interplane 
links, and connect two spinon sites i and j along various 
paths in the planes with these ladders; the fictitious flux 
connecting the spinous will then destroy the correlations 
across all the rungs on each ladder. 

The simplest estimate for the strength of the confin- 
ing force is to assume a circular droplet of decorrelated 
links between the spinous, with a diameter given by the 
spinon-spinon distance r. (This is presumably not a valid 
approximation for large spinon separations, but may be 
reasonable for the ground state of the oscillator.) It yields 
a harmonic potential 
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The ground state energy of the spinon-spinon bound 
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FIG. 5. Radial wave functions for the lowest energy eigen- 
states of the linear oscillator in two dimensions. 



state is hence given by 
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where we have used §) and 1(8*8^)^1 w 0.62 J^/ J. This 
estimate for the spin gap [IJ] in YBa2Cu306+a; is di- 
rectly related to the optical magnon gap in the ordered 
antiferromagnet YBa2Cu3 06.2, which has been measured 
by inelastic neutron scattering |15|: 
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which yields a spin gap of 27 mcV. The charge gap is 
just the gap to create a spinon-holon bound state: sub- 
stituting J = 120 meV, J± = 10 meV t = 500 meV, and 
tcS — 0.9 1 for YBa2Cu306+2: into (^), we obtain 31 meV. 

The microscopic details of the chirality liquid ||l^ yield 
a linear potential as an estimate for the confining force, 
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where Eqq is the spin gap given by the ground state en- 
ergy of the string oscillator (^ with string tension F2d in 
two dimensions. This oscillator has to be solved numeri- 
cally; the solutions are 
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with rp = l/(2mi.cd^2d)^^'^, A„f as listed in table |, and 
(j)ne{r) as plotted in fig. ||. 

Solving (P) and (|l|) for £^00 yields 
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where we have used and |(SiSj)_L| ~ 0.62J_l/J once 
more. This yields a spin and charge gap of 25 meV and 
28 meV for YBa2Cu306+x. 

It should of course be born in mind that these num- 
bers are only rough estimates; many important details, 



including the d-wave symmetry of the superconducting 
order parameter, have not been taken into account here. 

5. I wish to conclude with a few general observations 
relating the circle of ideas drawn above to high-Tc super- 
conductivity. 

(a) The fictitious flux confinement forces are not re- 
sponsible for the pairing of the charge carriers in a con- 
ventional sense, as the charge superfluid forms an integral 
part of the spin liquid within the planes; in particular, 
the interplane correlations are not destroyed by the pres- 
ence of charge carriers in the superconducting state at 
T = Q. The charge gap estimated above, however, cor- 
responds roughly |l6[ to the gap required to create an 
"unpaired" hole, which does not participate in the su- 
perfluid. The confinement forces are thus responsible for 
the robustness of the pairing, including its persistence at 
finite temperatures, but not for the pairing itself |l7|| . 

(b) The strength of the magnetic pairing force is deter- 
mined by the the magnetic energy stored in the structure 
between the layers, or more precisely by the amount this 
energy is raised by the fictitious flux, which is only in 
multilayer materials directly related to the coupling be- 
tween the layers. The mechanism hence allows for strong 
pairing forces in single layer materials; we expect those, 
however, to be crucially dependent on the magnetic prop- 
erties of the structure between the planes. 

(c) My thinking on the spin gap phase coincides with 
that of Emery and Kivelson |18|. 

(d) The superconductivity in Cgo is, according to my 
thinking, due to Dirac quantization of the fictitious flux 
through the surface of each molecule ]l9| , p0| , p^ . 
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minating discussions. This work was supported through 
NSF grant No. DMR-95-21888. Additional support was 
provided by the NSF MERSEC Program through the 
Center for Materials Research at Stanford University. 
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